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ABSTRACT 

Due to its proximity to the Earth and its nearly edge-on geometry, the Orion Bar pro- 
vides an excellent testbed for detailed models of the structure of Hn regions and the 
surrounding photon-dominated regions. In the present study, a self-consistent model 
of the structure of the Orion Nebula in the vicinity of the Bar is built under the as- 
sumption of approximate ionization, thermal, and hydrostatic equilibrium. It is found 
that a fairly simple geometry is able to describe the surface brightness profiles of the 
emission lines tracing the ionized Hn region with a remarkable accuracy, independent 
of the prescription adopted to set the magnetic field or the population of cosmic rays. 
Although we consider different scenarios for these non-thermal components, none of 
the models is able to provide a fully satisfactory match to the observational data for 
the atomic layer, and the predicted column densities of several molecular species are 
always well above the measured abundances. Contrary to previous studies, we conclude 
that a more elaborate model is required in order to match all the available data. 

Key words: ISM: HII regions - ISM: PDR - ISM: individual (Orion Nebula, M42, 
NGC 1976) - ISM: individual (Orion Bar) 



1 INTRODUCTION 

Hn regions are extended, low-surface brightness, diffuse neb- 
ulae of photoionized gas. They are associated with regions of 
ongoing star formation, most commonly found in the disks of 
spiral galaxies, where young and bright stars provide the co- 
pious amount of ionizing ultraviolet (UV) radiation required 
for such regions to exist. The spectrum of these nebulae is 
mainly composed of hydrogen recombination lines and for- 
bidden lines of ions of common elements, superimposed on 
a weak continuum. Hn regions can be used to trace star 
formation from the solar neighborhood to the high-redshift 
Universe, and their spectrum provides invaluable informa- 
tion about the ionizing population of massive stars, as well 
as the physical conditions of the interstellar medium. 

The central source of such a region can be one or sev- 
eral Population I stars of type O or early B, with effective 
temperatures between 3 and 5 x 10 4 K, that emit a large 
number of photons with energies higher than the ionization 
potentials of hydrogen and helium. Although these two ele- 
ments are by far the most abundant ones, metal lines play 
an important role because they provide the principal cooling 
mechanism. At any point in the nebula, the degree of ion- 
ization is determined by the equilibrium between electron 
capture and photoionization. Ejected photoelectrons carry 
the excess energy of the photon as kinetic energy, and they 
contribute through electron-electron and electron-ion colli- 



sions to maintaining a Maxwellian velocity distribution with 
typical electron temperatures between 5000 and 20000 K. 

The limit of the Hn region, an ionization front that 
tends to expand into the surrounding neutral gas at sub- 
sonic velocities, is usually approximated as a Stromgren's 
sphere with typical radii of the order of parsec. The Hn re- 
gion is surrounded by a neutral, photon-dominated region 
(PDR), where the stellar UV radiation heats the gas and 
partly dissociates the molecular hydrogen (hence the name 
photo-dissociation region, which is also often used). Pho- 
tons escaping far into the outer molecular region are also 
absorbed by the dust present in the interstellar medium, 
which gets thus heated to about 100 K, and re-emitted as 
an infrared continuum. 

Notwithstanding with this simple picture, observations 
of many Hn regions show signatures of dense neutral conden- 
sations scattered in the ionized zone and turbulent motions 
within the gas with velocities of the order of 10 km s _1 . 
The geometry of these nebulae is anything but spherical, 
and theoretical models of their internal structure ought to 
be constructed on case-by-case basis. 

In this work, we attempt to build a self-consistent model 
of the Orion Bar. The Orion Nebula (M42, NGC 1976) is 
part of th e Orion Molecular Cloud Complex, at a distance 
of 437 pc (jHirota et alj|2007l ). and it has been studied exten- 
sively over the years across the entire electromagnetic spec- 
trum. Most of the ionizing radiation comes from the star 
8 1 Orionis C, located roughly at the centre of the region. 
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The Orion Bar, situated at approximately 0.235 pc from the 
central star towards the south-east, is a dramatic example of 
the interface between the ionized gas and the PDR. Due to 
its convenient orientation, allowing an almost edge-on per- 
spective, and the rich variety of observations available in the 
literature, the Orion Bar constitutes an ideal laboratory for 
testing the PDR phy sics. 

Recent studies db'Dell et al.l 120091 : IShaw et al.1 120091 : 
IPellegrini et al.l 120091) highlight the importance of cosmic 
rays and magnetic fields in determining t he structure of 
the P DR. In particular, it has been argued (|Pellegrini et al.l 
l2009h that a relatively strong magnetic field, as well as a pop- 
ulation of cosmic rays in equipartition (and thus a density 
that is much higher than the average galactic background), 
must be present in the Bar in order to reproduce the ob- 
served surface brightness profiles of the H2 line at 2.121 fim 
and the 12 CO(J =1 — 0) emission. 

The presence of an important magnetic field is sup- 
ported by polari zation observations of the Orion Nebula 
l|Schleuning||l998T l and the detection of a magnetic field in 
Orion's Veil with an intensity that is at least an order of 
magnitude higher than the typical value s measured in th e 
cold neutral medium of the Milky Way (|Abel et al.ll2006l ). 
but the enhancement of the cosmic ray contribution above 
the Galactic background is much more poorly constrained 
from the observational point of view. A high cosmic-ray den- 
sity is simply introduced as an additional heating source that 
can act much deeper into the molecular cloud than the pho- 
tons from the central star. As a result, the temperature of 
the outer regions is considerably higher, providing a better 
fit to the observed surface brightness profiles. 

The present work represents an additional step towards 
a self-consistent model of the internal structure of the Orion 
Bar. As previous studies, it is based on the assumption that 
the region is in approximate hydrostatic equilibrium: the 
outward acceleration driven by the radiation field and the 
pressure of the hot, ionized gas creates an expanding wind 
that compresses the surrounding medium - and the mag- 
netic field coupled to it - until the magnetic pressure is able 
to halt the process. Radiative transfer across the different 
gas phases is solved by the photo ionization code Cloudy (last 
described in lFerland et al.lll99&t ). including all the relevant 
processes affecting atoms, molecules, and dust grains, and 
the predicted line intensities are computed by integrating 
the physical properties of the gas along the li ne of sight, fol- 
lowing an approach very similar in spirit to IMorisset et all 
l|2005i ). 

The main improvements with respect to previous work 
are the inclusion of a detailed treatment of the gravita- 
tional acceleration and a more elaborate description of the 
three-dimensional geom etry of the system. As shown in 
lAscasibar fc Di'azl (|2010l ), the gravitational force (both the 
mass of the central object as well as the self-gravity of the 
gas) plays an important role in the outer regions, setting the 
total extent of the molecular layer and eliminating one free 
parameter of the model. In addition, we propose a simple 
parameterization of the geometry of the Hn region that is 
able to provide a reasonable fit to the emission line data. 
Emission from the atomic layer, though, as well as the col- 
umn densities of several molecular species, remain difficult 
to reproduce for any model, with or without gravity. 

The details of our photoionization models are thor- 



oughly described in Section[2] and their predictions are com- 
pared with observational data in Section [3] Section [4] is de- 
voted to the physical interpretation of our results, and a 
brief summary and outlook are provided in Section [5] 



2 A SIMPLE MODEL OF THE ORION BAR 

Our model of the Orion Bar is based on the assumption that 
the gas is in ionization, thermal, and hydrostatic equilibrium 
at every point. Under these conditions, the physical prop- 
erties of a cloud with plane-parallel or spherical symmetry 
can be e fficiently computed with the plasma physics code 
CloudjQ (|Ferland et aI.|[T998r i. a spectral synthesis program 
designed for the study of low-density environments that are 
ionized by an external radiation field. In this section, we 
describe the parameters of our photoionization models and 
discuss how the condition of hydrostatic equilibrium and the 
non-spherical geometry of the Orion Nebula in the vicinity 
of the Bar may be implemented. 

2.1 Model parameters 

On input, Cloudy requires the user to specify the shape and 
the intensity of the incident radiation field, the gas density, 
its chemical composition, and the geometry of the cloud. For 
consis tency with previous studies, we follow IPellegrini et al.l 
(|2009T) and represent the incident con tinuum by the sum 
of the cosmic microwave background, a lKuruczl (|l979l ) stel- 
lar atmosphere with temperature T = 39600 K and ioniz- 
ing luminosity Q(H) = 10 49 photons s _1 , and a thermal 
bremsstrahlung component with temperature T = 10 6 K 
and luminosity L = 1 32 6 erg s _1 in the 0.5 — 8 keV band 
|Feigelson et al.ll2005l ). The gas density at the illuminated 
face is set to no = 10 3 ' 2 cm~ 3 , and a chemical com position 
appropriate for Orion Nebula (|Baldwin et al.|[l99ll ). includ- 
ing dust grains and polycyclic aromatic hydrocarbons, is 
used. 

Unless otherwise specified, Cloudy assumes that the gas 
density stays constant over the cloud, but many density and 
pressure laws can be used instead. We model the Orion Bar 
by enforcing hydrostatic equilibrium with the constant pres- 
sure command. Contrary to what the name suggests, this 
option does not keep the total pressure fixed, but adjusts its 
value so that the total acceleration vanishes at every point, 

and therefore 

/• r 

P(r) = Po + / p{x) a(x) dx (2) 
Jo 

where Pq denotes the initial pressure at the illuminated 
face, and the acceleration to be balanced includes the terms 
due to the absorption of photons (pushing the gas away 
from the central star) and the gravitational force that 
pulls the whole cloud towards the centre. The latter is 
included by means of the gravity command, described in 

All our calculations have been preformed with version C08.00 
of the code. The current stable release is publicly available at 
http:/ /www. nublado.org 
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lAscasibar fc Dfaj (|20ld ). The gravitational acceleration is 
calculated as g(r) = — A-kGM(t) /r 2 , assuming spherical 
symmetry, and the mass inside radius r includes the con- 
tributions of both the gas and the stars in the cloud. The 
distribution of the gas mass is computed self-consistently by 
the code, M gas (r) = p(x)4irx 2 dx, and the 1 Ori C multi- 
ple sy stem is modeled as a point mass of 50 Mq (|Kraus et al.l 
12007 ) located at the centre of the region. 
The total pressure 



P(r) 



+ Pur 



Plin 



(3) 



that appears in equation ([2]) is the sum of the thermal 



pressure of the gas P g; 



nkT, the terms P ra m = P^win 



and Pturb = P^turb/S induced by the uniform and turbu- 
lent motions, respectively, the ma gnetic pressure Pnag = 
B 2 /8n, and the contribu tion P incs (|Ferland fc Elitzur|[l98i ; 
lElitzur fc Ferlandlfl98r3 ) of the trapped emission lines. Al- 
though winds are not considered in our model, it includes 
a small turbulent velocity field of 2 km s _1 . Regard- 
ing the magnet i c field , we consider the same scenarios as 
IPellegrini et all |2009t) . In the gas pressure model, only the 
thermal pressure of the gas and the turbulent pressure terms 
are taken into account. The magnetic pressure model adds 
a tangled magnetic field whose intensity at the illuminated 
face of the cloud is Bo = 8 p,G, and its equation of state 
is given by 7 = 2 (i.e. B/Bo = n/no). In both scenarios, 
the density of cosmic rays is constant and equal to the aver- 
age galactic background. The enhanced cosmic rays model 
assumes the cosmic-ray density to be in equipartition with 
the magnetic field, resulting in a much higher abundance of 
relativistic particles. 

For each scenario, we ran a grid of Cloudy models with 
initial distances from 9 1 Ori C ranging from 0.05763 to 
0.5763 pc in logarithmic steps of 0.01. A sample Cloudy 
script used in the preparation of our model grids is shown 
in Appendix fS] Each model saves the emissivities of two 
emission lines from the ionized region ([Sll] A A 6716+6731 
and Ha) and one emission line produced in the atomic layer 
(H2 2.121 /im), the extinction towards the end of the cloud, 
and the volume densities of six molecular species (CO + , SO, 
CN, CS, SO + , and SiO), all of them as a function of depth 
into the cloud (i.e. the distance from the illuminated face). 

2.2 Geometry 

One of the aims of the present work is to show that the over- 
all geometry of the Orion Nebula plays an important role on 
its physical conditions and observable properties, which can 
be exquisitely probed in the region near the Bar due to its 
privileged orientation. Previous studies have constrained the 
geometry of the cloud using the surface brightness profile of 
the emission lines associated to the Hll region. More specif- 
ically, the [Sll] A A 6716+6731 line displays a sharp peak at 
the interface between the ionized and the neutral layers, thus 
providing an e xcellent tracer of the position of the ioniza tion 
front (see e.g. iBaldwin" et al.lll99ll ; fWen fc Q'delllll995h . 

In the vicinity of the Bar, the projected distance be- 
tween the maximum of the observed emission and the central 
star is 111 arcsec (0.235 pc, assuming a distance of 437 pc). 
Approximating the Bar as a plane-parallel slab of thickness 
h, located at a distance Po from d 1 Ori C with an incli- 
nation angle /3 with respect to the line of sight, the data 




Figure 1. Geometry of the illuminated face of the cloud, de- 
scribed by the parameters Rq, a, (3, 8, and h (see text for a de- 
tailed explanation). The x- and y-axes represent the offset with 
respect to 8 1 Ori C and the coordinate along the line of sight, 
respectively. 



suggest the values h ~ 0.115 pc, Po ~ 0.114 pc, and j3 ~ 7° 
^Pellegrini et alj|2009h . 

We consider a very similar layout, depicted in Figure [1] 
that also includes the contribution of the adjacent regions of 
the Orion Nebula. The interior region, towards the central 
star, is approximated as a sphere of radius Ro centered at 
the location of 1 Ori C, which we take as the origin of 
coordinates. The transition between the spherical region and 
the Orion Bar takes place at an angle a with respect to 
the plane of the sky, which is an additional free parameter 
of our model. The Bar itself, as well as the outer region, 
are modeled as straight lines whose angles with respect to 
the line of sight are denoted by the parameters /3 and S, 
respectively. The length of the Bar region is specified by 
the parameter h, whereas the outer region is assumed to 
continue well beyond the observed field. 

The region interior to this curve is empty, or, more pre- 
cisely, fi lled by a hot, tenu ous gas that can be detected in 
X-rays (|Giidel et al.l I2008T 1 and provides thermal pressure 
support for the warm (T ^ 10 4 K) dense gas that is respon- 
sible for the optical and infrared emission in which we are 
interested. In Cloudy parlance, the curve depicted in Fig- 
ure[T]defines the illuminated face of the cloud. For the outer 
region to receive direct illumination from 8 1 Ori C, as as- 
sumed by our models, it must not be behind the shadow 
cast by the Bar (angle ip in Figure [1]). Thus, we chose the 
parameterization 



(4) 



where < ?y < 1. 

This geometry is implemented as a summation over a 
grid of Cloudy models with spherical symmetry and different 
distances between the illuminated face of the cloud and the 
central star (light dashed lines in Figure [1] illustrate the sep- 
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aration between models). The spherical region corresponds 
to a single Cloudy model, fully characterized by the value of 
Ro, while both the Bar and the outer region involve many 
models each. 

The column density of a given species i is given by the 
integral along the line of sight 



Ni 



m(x,y) dy 



where the volume density m(x, y) is evaluated from the out- 
put of the Cloudy model that is appropriate for each posi- 
tion. Surface brightness profiles of the emission lines have 
been computed as 



Si{x) 



£i( x i y) ^Q-0AAi(x,y) 



4tt 



dy 



(6) 



where 6i(x, y) denotes the emissivity per unit volume of the 
line, and the amount of dust extinction 



Ai(x,y) = 



dy 



j{x,y) dy 



includes all the foreground material between any given point 
and the observer. The differential extinction ^f-(x,y') is 
estimated from the radial increment of the total extinction 
in the V band output by Cloudy and then converted to the 
rest wavelength of the line. 
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Figure 2. Observed surface brightness profiles of the 
[Sn] A A 6716+6731 and Ho emission lines (dots), compared to the 
theoretical model predictions. The gas pressure, magnetic pres- 
sure, and enhanced cosmic rays models have been plotted as dot- 
ted, dashed, and solid lines, respectively, but they lie virtually on 
top of each other. 



3 COMPARISON WITH OBSERVATIONS 
3.1 Hn region 

As will be shown below, the observational properties of the 
ionized region (more precisely, the surface brightness profiles 
of its emission lines) are not very sensitive to the details of 
the magnetic field or the prescription adopted to establish 
the population of cosmic rays. The exact position of the 
ionization front depends on a combination of the gas density 
at the illuminated face, the intensity of the ionizing radiation 
from 6 Ori C, and the distance to the Orion Nebula. Once 
these (degenerate) parameters are specified, the shape of the 
emission profiles is entirely determined by the geometrical 
configuration of the system. 

In our model, this geometry (the distance from the cen- 
tral star to the illuminated face) is described by the values 
of five free parameters: the radius of the inner region Ro, 
the transition to the Bar a, its inclination with respect to 
the line of sight /3, its length h, and the angle 5 - or, equiv- 
alently, the parameter r\ defined in expression (J3J> _ that 
defines the orientation of the outer region with respect to 
the line of sight. We estimate the values of these parame- 
ters by fitting the observed surface brightness profiles of the 
[Sn] A A 6716+6731 A and Ha emission lines. 

The surface brightness profile of the [Sn] A A 6716+6731 
line across the ba r has been obtained fro m two narrow-band 
images, taken by iPellegrini et aL | (|2009h with the Southern 
Astrophysical Research Telescope, centered at A = 6723 and 
6850 A with bandwidths of 45 and 9 5 A, respectively. Fo r 
the Ha line, we use the observations of lWen fe O'deil l|l995h . 
Both data sets were constructed as continuum-subtracted 
averages over 20 arcsec-wide swathes, using similar cuts. 



Since these data have already been corrected for dust ex- 
tinction, we set Ai = when modelling both emission lines. 

The best-fitting values of Ro, a, /3, r\, and h have 
been found by m eans of the FiEstAS sampling technique 
I Ascasibarl |2008| ) . a Monte Carlo integration scheme based 
on the Field Estimator f or Arbitrary Spa ces (FiEstAS; 
lAscasibar fc Binnevl 120051 ; lAscasibarl l2010h . In order to 
quantify the quality of the fit to the observational data, we 
compute the reduced x 2 as 

2 _ 1 [S ohs (xj) - S mo dclfe)] 2 , , 



where S h a {xj) denotes the observations at a projected dis- 
tance Xj from the central star, S mo d e i (x 3 •) are the cor- 
responding model predictions, and the sum over the in- 
dex j corresponds to the iV bs observational data points, 
whose errors have been characterized by a standard devia- 
tion a 3 ■ — 5 x 10~ 14 and 5 x 10~ 13 erg s _1 cm -2 arcsec -2 
for the [Sn] A A 6716+6731 and the Ha lines, respectively. 

The results are shown in Figure where the surface 
brightness profiles obtained for Ro = 0.122 pc, a = 28.3°, 
f) — 2.3°, h = 0.05 pc, and 8 — 45.15° are compared with 
the observational data points. Our three scenarios for the 
magnetic field and cosmic rays yield exactly the same pre- 
diction for the emissivities within the ionized region of the 
nebula. Although the value of the best-fit \ 2 — 0-62 is some- 
what arbitrary, since it depends on the adopted ai, we judge 
from Figure [2] that the proposed geometry is able to provide 
a reasonable description of the Orion Bar up to the ioniza- 
tion front, regardless of the assumptions made concerning 
the non-thermal components. 
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3.2 Atomic region 

In order to test the ability of our models to describe the 
structure and physical properties of the atomic layer, we 
will now focus on the H2 S(l — 0) transition at 2.121 /xm, 
comparing the model predictions with t he observational 
data reported in Ivan der Werf et~ai] (| 19961 ). The spatial cut 
from which the H2 data were obtained is different from 
that of the [Sn] and Ha lines. Although the [Sill pro- 
file shows little variation (see e.g. iHennev et all 20051). it 
has been shown (Ivan der Werf et al.l 19961; lYoung Owl et al.l 



l200d ; IWalmslev et al.ll2000l ; lAllers et alj|2005l ) that the sur- 
face brightness of the H2 line may change considerably when 
using different cuts perpendicular to the Bar, which intro- 
duces some uncertainty in the comparison. Moreover, the 
data represent a summation parallel to the bar, whereas the 
H2 emission is actually concentrated in rather thin filaments 
with typical widths of ~ 8 arcsec, about 0.016 pc. 

Another potential issue is the estimation of the amount 
of dust extinction in the models. On the one hand, our 
models may overestimate the value of Ay, per haps due to 
the a ssumed cross-section for dust attenuation |AUers et al.l 
120051 ) and/or an excessive amount of dust in the atomic 
region. On the other hand, the conversion between the visual 
extinction Av returned by Cloudy and the extinction at the 
infrared wavelength A = 2.121 fim of the H2 emission line 



A 



2.121 fim 



0.14 A 



v 



(0) 



may be obtained by applying a ICardelli et ail (| 19891 ) red- 
dening curve with R — 5.5, in accordance with observations 
of the Orion Nebula. This estimate is based on the observed 
reddening curve of stellar spectra, where the extinction in- 
cludes both the absorption and scattering of light by the 
dust grains. For an extended source, like the Orion Nebula, 
a large fraction of the photons (those that are scattered a 
small angle away from the light of sight) will be compensated 
by similar small-angle events affecting nearby rays, and the 
effective scattering opacity 



(1-5) 



(10) 



will be given by the product of the total scattering cross- 
section a s and the grain asymmetry factor g, defined as the 
average {cos 8) over the angle 9 between the incident and 
the scattered photon. Since the wavelength dependence of 
the asymmetry factor is determined by the chemical compo- 
sition, shape, and size distribution of the dust grains, and 
these variables may vary with the spatial location within the 
nebula (they depend, for instance, on the gas density and 
the high-energy radiation field), there is some uncertainty 
associated to equation ((9}. Nevertheless, an independent 
estimation of the nebular extinction in the Orion Nebula 
based on the hydrogen recomb ination lines from the Balmer, 
Paschen, and Brackett series l|Bautista et al.ll"l995h is com- 
patible with A2. 121pm = 0.14 71V, suggesting that the effects 
of scattering are not very important in this case. 

The predicted surface brightness profiles of the 
H2 2.121 fj,m line for the gas pressure, magnetic pressure, 
and enhanced cosmic rays models are shown on the different 
panels in Figure [3] It is worth noting that, while the emis- 
sion lines of the ionized region have been used to constrain 
the free parameters that set the geometry of the illuminated 




x Lpcj 

Figure 3. Surface brightness profiles of the H2 2.121 (im emission 
line predicted by the gas pressure (top) . magnetic pressure (mid- 
dle), and enhanced cosmic rays (bottom) models. Dotted lines 
show the total, uncxtinctcd emission. Dashed and solid lines cor- 
respond to A 2 .i2i M m = 0.028 Ay and A3.121.pm = 0.14 Ay, re- 
spectively. 



face, the surface brightness profile of the H2 line is a genuine 
prediction of the models. 

Our results suggest that none of the scenarios we have 
considered is able to reproduce the observed emission at the 
quantitative level. In particular, all models underestimate 
the maximum intensity of the H2 line by more than a fac- 
tor of three. We have verified that the discrepancy between 
models and data is not a consequence of the prescription 
adopted for dust extinction. In the gas pressure scenario, 
the separation between the peaks of the [Sn] and the H2 
emission is always too small, whereas the magnetic pressure 
model fails to explain the observed surface brightness even 
for A2. 121 f*m = 0. The enhanced cosmic rays model seems 
slightly more promising, but the precise amount of dust ex- 
tinction is critical. In addition, the density gradient in the 
molecular region should be much steeper than predicted by 
our gravitational models in order to reproduce the observed 
lack of emission outside the main peak. 



3.3 Molecular region 

The abundances of different molecular species pro- 
vide a direct probe of the extent and the physi- 
cal conditions of t he outer, molecular layer. As in 
IPellegrini et~ai] (|200Sh . we compare the column density pro- 
files predicted by the models wit hobser vat i ons of CO + 
JStoerzer et al.l 1 19951; IFuente et alJl2003D. CN (ISimon et all 



Il997f). SO+ llFuente et al.ll2003l). SO jjansen et all Il995l ). 



CS Isimon et alj Il997l ; iHogerheiide et alj Il995l )~and SiO 
(|Schilke et alj|200ll ). 

Observed column densities should, however, be treated 
with some caution, given the ma ny assumptio n s invo lved 
in their derivation. For instance, IFuente et al.l (|2003l ) as- 
sume a constant rotational temperature of 10 K for all the 
molecular levels, and all lines are considered to be opti- 
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Figu re 4. Observed column den sities of CO+ jStoerzer et al]|l995t IFuente et al.|[20o3), CN llSimon et al]|l997h . SO+ llFuente et al.l 
l2003t ). SO {j arise n et al Jl995h . CS {| imon et al Jl997tlHogerheiide et aljl995h . and SiO jSchilke et alj200ll ). compared to the predictions 
of the gas pressure (dotted), magnetic pressure (dashed), and enhanced cosmic rays (solid lines) models. 



cally t hin. For the CN and CS column densities. ISimon et al.l 
combined their ob servations with prev ious data from 
IFuente et all \ 19961 ) and Ivan der Werf et al.l (|l996h . They 
assumed kinetic temperatures of 40 — 100 K and used CS 
rates for the unknown CN collisional rate coefficients. CN 
optical depths were estimated by an escape probability code, 
assumi ng a single compon e nt. Ho wever, the CS column den- 
sity of iHogerheiide et ah! (|l995l ) was derived using statis- 
tical equilibrium calculations that assumed dense clumps 
to be embedded in a lower density medium, and the sa me 
type of analysis was employed by Ijansen et al. | (|l995h to 
in fer the column dens ity of SO. The SiO column densities 
of ISchilke et~ai1 |200ll ) were obtained from observations of 
three transitions, assuming th e same temperatu r e in the 
molecular region, 85 ± 30 K, as IHogerheiide et al.l (|l995l ). 

For this reason, as well as several factors that may affect 
the accuracy of the models (see the discussion below), the 
comparison between model predictions and observational 
data is far from straightforward. With this caveat in mind, 
our results for the six molecular species considered are shown 
in Figure U 

The CO + column density, shown on the upper left 
panel, is reasonably well reproduced by the enhanced cos- 
mic rays model. The intensity of the peak, as well as the 
offset at which it occurs, are consistent within the uncer- 
tainties, and the model also explains the extended wings of 
the column density profile. However, the gas pressure and 
magnetic pressure models predict about a factor of 100 
smaller colum n densities than observ ed, in agreement with 
the findings of iPellegrini et ah! (|2009T ). On the contrary, for 
the CN molecule, represented in the top middle panel, these 
models are r emark ably close to the observational data of 
ISimon et alj (fl997h . while the enhanced cosmic rays predic- 
tions are now two orders of magnitude above the observed 
values. 



In both cases, the shape of the profiles is arguably a 
little too flat. Such a slow increase of the column density 
when moving away from the central star into the molecular 
region is even more evident for SO + (shown on the top right 
panel), whose abundance is bracketed by the predictions of 
the gas/magnetic pressure and enhanced cosmic rays mod- 
els, and the SO molecule, plotted on the bottom left panel, 
where the models with Galactic cosmic rays overestimate the 
observed column density, and none of them reproduces the 
observed gradient. The same occurs for CS on the bottom 
middle panel, where all models are more than two orders 
of magnitude above the data (in contrast with the findings 
of IPellegrini et all [2009. where the magnetic model gives a 
good match), and SiO, on the bottom right panel, where the 
discrepancy reaches four orders of magnitude. 



4 DISCUSSION 

To summarize, our models of the Orion Bar seem to pro- 
vide an excellent fit to the Hn region, but they do a much 
poorer job for the atomic and molecular layers. On the one 
hand, they tend to overestimate the molecular column den- 
sities under all the scenarios considered for the intensity of 
the magnetic field and the population of cosmic rays. On 
the other hand, the comparison between the predicted and 
observed intensities of the H2 emission line hints that the 
models also overestimate the amount of dust in the atomic 
layer. 

One possible reason is, of course, that our models are 
too simple to describe all the relevant features of the Orion 
Nebula. Given the computational cost of building one of our 
model grids, a full exploration of the parameter space (where 
the cost increases exponentially with the number of free pa- 
rameters) is well beyond the scope of the present work. Nev- 
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Figure 5. Number density of ionized, atomic and molecular 
hydrogen (shown as blue, black, and red lines on the top panel, 
respectively), as well as electron temperature (bottom panel) for 
the enhanced cosmic rays case (solid lines) and a model with 
Bo = 50 fiG and 7 = 1 (dashed lines). 

ertheless, we have run several single Cloudy models to test 
whether a better agreement with the observations could be 
achieved by varying the density of the gas at the illuminated 
face, the intensity and equation of state of the magnetic field, 
or the amount of mass in stars. Alternatively, it is also pos- 
sible that the proposed models are incomplete, not because 
they lack the required flexibility, but because they do not 
include all the relevant physics. We briefly discuss the va- 
lidity of some of the main assumptions of the model - most 
notably, the approximation of hydrostatic equilibrium. 



of the atomic phase may be adjusted by varying Bo or 7, 
with larger values of either quantity yielding smaller den- 
sities and more extended regions for both the atomic and 
molecular layers. 

One could try to obtain a good fit to the observational 
data by modifying Bo and 7 at the same time. A rigor- 
ous analysis would require the computation of an extensive 
model grid, but we have experimented with several combi- 
nations, and none of them seemed to be particularly promis- 
ing. As an example, we show in Figure the density and 
temperature structure of a model with Bo = 50 [iG and 
7 = 1, consistent with observations of molecular clouds 
( Hci les fc Crotchet! [20051 ) and radiative MHD simulations 
of th e expansion of magnetised Hn regions (|Arthur et alj 
l201lf ). The predictions of this model are almost indistin- 
guishable from the standard enhanced cosmic rays scenario. 
More generally, all parameter choices where the density was 
high enough to produce significant H2 emission at 2.121 /im 
led also to a very high dust extinction and thus to surface 
brightness profiles similar to those shown in Figure [3] 

However, these tests do not exhaust the available de- 
grees of freedom. In particular, there is no reason why the 
gas density and the magnetic field should be uniform ac ross 
the illuminated face. According to IWen fc O'delil (|l995l . see 
their Figure 4) , the electron density near the ionization front 
ranges from ~ 2000 to ~ 7000 cm -3 along a cut perpen- 
dicular to the Bar, and observations frequently show mul- 
tiple peaks in the H2 emission as one moves away from 
the Trapezium across the Bar. In many cases, some of 
these peaks occur inside the main ionization front (closer 
to the Trapezium tha n the [Sll] peak, see e.g. Figure 4 of 
IWalmslev et~aLll2000l ). providing further evidence that the 
Bar geometry is quite complex, with multiple overlapping 
filaments. Any realistic model of the Orion Nebula should 
account for these changes in the gas density encountered by 
different rays, which are probably associated with similar 
variations in the magnetic field. 



4.1 Gas density and magnetic field 

To a certain extent, the initial gas density of the models 
can be thought of as a global scaling factor. Roughly speak- 
ing, column densities scale approximately proportional to 
no, emissivities increase as n 2 ,, and distances are propor- 
tional to rig 1 . Although there is some freedom in the choice 
of the optimal value of this parameter, especially if one al- 
lows for errors in the estimation of the distance to the Orion 
Nebula, the observed line ratios of several emission lines im- 
pose severe constraints on the ionized gas density. Moreover, 
once the geometry is adjusted to fit the surface brightness 
profiles in the Hn region, the predictions for the other quan- 
tities are remarkably insensitive to the precise value of no. 

As explained in Section 12.11 the magnetic field is de- 
scribed by the equation of state 




with Bo — 8 /iG, no = 10 3 ' 2 cm~ 3 , and 7 = 2 in all our 
models. In general terms, the intensity of the magnetic field 
regulates the density contrast between the different phases. 
It has a negligible effect on the ionized region, whose density 
is basically set by the value of no . The density and thickness 



4.2 Stellar mass distribution 

In addition to the gas density at the illuminated face and 
the equation of state of the magnetic field, we have also 
considered different prescriptions to account for the gravita- 
tional acceleration associated to the stellar component. The 
influence of other stars in the Trapezium has been studied 
by changing the central mass from 50 to 200 Mq (see e.g. 
iHerbig fc Terndrup|l98^ ; lHuff fc Stahleil200d ). and we have 
also considered a continuous mass distribution of the form 

M(r)=M«— (12) 

with M* = 2000 M and r* = 2 pc to repre- 
sent the whole stellar populat ion of the Orion Complex 
(|Hillenbrand fc Hartmaimlll998! '). 

The effect of the gravitational acceleration is shown in 
Figure [6] where we compare the gas density, temperature, 
and pressure profiles predicted by the enhanced cosmic rays 
model in the absence of gravity, our standard implemen- 
tation (self-gravity of the gas plus d 1 Ori C), and a model 
where the stellar mass is represent ed by equation (|12ll . As i n 
the plane-parallel case discussed in lAscasibar fc Dfazl (|201dh . 
the physical conditions of the gas in the ionized and neutral 
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Figure 6. Results obtained in the absence of gravity (dotted 
lines), considering 6 1 Ori C as a point mass of 50 Mq (solid 
lines), and using equation l i 1 21 ) to model the total stellar mass 
distribution (dashed lines). Top panel: number density of ion- 
ized (blue), atomic (black), and molecular (red) hydrogen. Mid- 
dle panel: electron temperature. Bottom panel: thermal (blue), 
magnetic (black), and turbulent (red) pressure components. 



layers are barely affected, fn the outer parts, both density 
and pressure decrease with radius in order to balance grav- 
ity, and the net effect is to effectively limit the extension 
of the molecular cloud. Although there is a qualitative dif- 
ference with respect to the non-gravitational case, in which 
the molecular layer is infinite, the prescription adopted to 
model the stellar distribution does not alter the results sig- 
nificantly. 

The importance of the gravitational term in equa- 
tion ([2} can be assessed by comparing the total gas pres- 
sure with the value obtained in the absence of gravity. From 
Figure [51 one sees that gravity plays a crucial role beyond 
~ 0.4 pc from the central star. At those distances, the mag- 
netic pressure dominates the total budget, but, due to the 
equation of state with 7 = 2, the turbulent component even- 
tually becomes dominant at larger radii as the gas density 
decrease^]. 

The escape velocity of the system, assuming the mass is 
distributed according to expression (|12[) . is about 2 km s _1 , 
comparable to the Alfven speed predicted by the models, the 
turbulent velocity dispersion that we assumed for the gas, 
and the observed velocity dispersion of t he molecular gas 
and the stars in the Orion Nebula Cluster (|jones fc Walked 



2 The ratio between magnetic and turbulent pressure would stay 
constant for 7 = 1. 



ll98Sl ; lTobin et al.ll2009l ). The predicted extent of the molec- 
ular cloud, of the order of 1 pc radius, is also consistent with 
the observed size. However, large-scale organised bulk mo- 
tions (champagne flow) and smaller scale chaotic motions, 
both at approximately sonic speeds (~ 10 km s~ ), have 
been observed in the ionized gas (see e.g. iGarcfa-Di'az et al] 
120081 ), whereas highly supersonic motions, of t he order of a 
few k m s _1 , are seen in the molecular gas (e.g. iTobin et al] 
120091 ) . suggesting that the Orion nebula is not a quiescent 
but a dynamically young object. 

4.3 Model assumptions 

First and foremost, the structure of the cloud may be com- 
pletely different if the approximation of hydrostatic equi- 
librium is relaxed. Indeed, the principal deficiency of our 
model is that it does not take into account gas motions 
and non-steady processes, which would change the emis- 
sivity of particular lines through the effects of advection 
and/or non-equilibrium chemistry (iBertoldi fc Draine|[T996l ; 
IStoerzer fc Hollenbachlll998l ; iHennev et al.ll2005l ) and mod- 
ify the pressure balance. The fact that the global model ge- 
ometry is neither spherically symmetric nor strictly plane- 
parallel means that there is also a formal inconsistency in 
the pressure balance, because it is only enforced radially but 
not in the transverse direction (between the different radial 
rays). However, this is probably a minor point in compar- 
ison with the assumption of a constant density across the 
illuminated face or the neglect of dynamics. 

In the restricted case of steady-state motions, the pres- 
sure balance along the ray could be dealt with by adding an 
extra term P ram = P^wind an< A adjusting the density so that 
the net acceleration does not vanish, but compensates for 
the chang es in radial velocity. Both observations an d simula- 
tions (e.g. iMellema et al]|2006l ; lArthur et ai]l201ll ) indicate 
velocities of the order of 10 km s _1 in the ionized region. 
Assuming a steady mass flux M implies the relation 

M 



(13) 



4-7rr 2 p 

which yields typical velocities ~ 1 km s _1 in the neutral lay- 
ers. The magnitude of the effect will thus be similar to that 
of changing the intensity of the magnetic field or its equa- 
tion of state. Although these arguments do not necessarily 
hold in the general, fully time-dependent case, we think it 
is unlikely that the presence of gas motions qualitatively 
changes the results obtained by imposing strict hydrostatic 
equilibrium. 

An independent check of the importance of non-steady 
effects is provided by the ratio of the sound-crossing time 
to the evolution timescales of the region of interest. Only if 
this ratio is small will a steady-state model be justified. The 
sound-crossing time 

Ar 

U = — (14) 

C s 

is of the order of 10 4 years in the ionized region (typical 
size Ar ~ 0.15 pc and sound speed c s ~ 15 km s _1 ) and 
the atomic layer (Ar ~ 0.015 pc, c s ~ 1.5 km s _1 ), but it 
is much longer in the outer, molecular layer (Ar > 0.2 pc, 
c s < 1 km s _1 ). Since these numbers are of the same order of 
magnitude as the evolution timescales of the corresponding 
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layers, estimated by dividing their hydrogen column densi- 
ties by the particle flux at the interfaces, the steady-state 
approximation seems to be only marginally valid. 

On the other hand, due to the relatively high ionization 
parameter, the gas cooling, hydrogen recombination, and 
H2 destruction timescales are much shorter than the sound- 
crossing time throughout most of the cloud, suggesting that 
the effects of non-equilibrium chemistry are very limited. 
The only place where the chemical and dynamical timescales 
may be comparable is the immediate vicinity (0.01 pc) of the 
dissociation front, which could affect the surface brightness 
profile of the H2 line at 2.121 /im. 

These simple calculations hint that the presence of gas 
motions would probably not modify the qualitative picture, 
but they must be taken into account in order to make a 
quantitative comparison with observational data. 

Regarding other model assumptions, the role of cosmic 
rays on the molecular chemistry is still unclear. Our results, 
most notably the abundances of CO + and SO + , hint that 
a large population of cosmic rays, well above the Galactic 
background but not as large as implied by equipartition ar- 
guments, is indeed present in the Bar, but the discrepancies 
found for the CN, CS, SO, and SiO molecules clearly in- 
dicate that cosmic rays alone cannot provide the ultimate 
solution to the problem. 

In addition, some atomic data (e.g. CS) are highly un- 
certain, and the theoretical predictions can not be more 
accurate than the atomic coefficients and reaction rates 
on which they are based. Another complication is that 
some of the observations (e.g. CS and SO) assume a two- 
component fluid, with dense clumps embedded in a more 
tenuous medium, rather than a single-phase gas. Further- 
more, Cloudy calculates the cloud structure by propagating 
all input continua from the illuminated face of the cloud 
towards the outer parts, but the gas in the molecular re- 
gion is also exposed to external radi ation sources, such as 
nearby stars (see e.g. the discussion in lO'Dell fc Harrisll2010l . 
concerning the role of 9 2 Ori A in the Orion Bar) or the ex- 
tragalactic UV background. 

In our view, the present discrepancies between the mod- 
els and the observational measurements highlight the com- 
plexity of Hn regions and, most notably, the surrounding 
PDR. Further theoretical effort is clearly required before we 
can claim a full understanding of the structure and chem- 
istry of these regions, as well as the main physical mecha- 
nisms that regulate them. 



All the photoionization models presented here have 
been computed with the plasma code Cloudy, under the as- 
sumption that the Orion Bar is in approximate hydrostatic 
equilibrium. In this scenario, the radius of the fully ionized 
region is determined by the gas density and the flux from 
the central star, whereas the thickness of the atomic layer 
is regulated by the intensity of the magnetic field. In the 
outer, molecular regions, the gravitational acceleration has 
to be balanced by a negative pressure gradient, and the total 
extent of the cloud is set by the action of gravity. 

Meaningful comparisons between model predictions and 
observational data require the assumption of an appropriate 
geometry for the system. In the particular case of the Orion 
Bar, previous studies have modeled the region as a plane- 
parallel slab of thickness h oriented with a small inclination 
angle /3 with respect to the line of sight. Including the con- 
tribution of the rest of the Orion Nebula has a significant 
impact on the surface brightness profiles of the emission lines 
produced in the ionized region, bringing them into excellent 
agreement with observations. 

However, none of the models presented here is able to 
reproduce all the observable properties of the Orion Bar. In 
particular, there is some mismatch in both the height and 
the width of the peak in the surface brightness profile of 
the H2 2.121 fim emission line from the atomic region, and 
a severe discrepancy in the molecular column densities. Al- 
though these results seem to be fairly robust with respect 
to variations in the model parameters, one possible solution 
would be to increase the complexity of the models - con- 
sidering, for instance, possible variations of the gas density 
across the illuminated face - and use the available data to 
constrain the additional degrees of freedom. 
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5 CONCLUSIONS 

Thanks to its proximity and orientation, the Orion Bar of- 
fers an exquisite view of the interface between the ionized, 
neutral, and molecular layers of the Orion Nebula. In the 
present work, we have developed self-consistent photoion- 
ization models of the Bar and compared their predictions 
with observational data available in the literature. In the 
light of previous works that have already constrained sev- 
eral physical ingredients that play an important role, such 
as the gas density, the magnetic field, or the abundance of 
cosmic rays, our study is devoted to explore the effects of 
the geometry of the region, including also the gravitational 
acceleration in the pressure balance of the gas. 



REFERENCES 

Abel N. P., Ferland G. J., O'Dell C. R., Shaw G., Troland 

T. H., 2006, ApJ, 644, 344 
Allers K. N., Jaffe D. T., Lacy J. H., Draine B. T., Richter 

M. J., 2005, ApJ, 630, 368 
Arthur S. J., Henney W. J., Mellema G., de Colle F., 

Vazquez-Semadeni E., 2011, MNRAS, pp 541-+ 
Ascasibar Y., 2008, Computer Physics Communications, 

179, 881 

Ascasibar Y., 2010, Computer Physics Communications, 
181, 1438 

Ascasibar Y., Binney J., 2005, MNRAS, 356, 872 
Ascasibar Y., Diaz A. I., 2010, MNRAS, 404, 275 



10 Y. Ascasibar et. al. 



Baldwin J. A., Ferland G. J., Martin P. G., Corbin M. R., 
Cota S. A., Peterson B. M., Slettebak A., 1991, ApJ, 374, 
580 

Bautista M. A., Pogge R. W., Depoy D. L., 1995, ApJ, 452, 
685 

Bertoldi F., Draine B. T., 1996, ApJ, 458, 222 
Cardelli J. A., Clayton G. C, Mathis J. S., 1989, ApJ, 345, 
245 

Elitzur M., Ferland G. J., 1986, ApJ, 305, 35 

Feigelson E. D., Getman K., Townsley L., Garmire G., 

Preibisch T., Grosso N., Montmerle T., Muench A., Mc- 

Caughrean M., 2005, ApJS, 160, 379 
Ferland G. J., Elitzur M., 1984, ApJ, 285, Lll 
Ferland G. J., Korista K. T., Verner D. A., Ferguson J. W., 

Kingdon J. B., Verner E. M., 1998, PASP, 110, 761 
Fuente A., Rodriguez- Franco A., Garcia- Burillo S., Martm- 

Pintado J., Black J. H., 2003, A&A, 406, 899 
Fuente A., Rodriguez-Franco A., Martin-Pintado J., 1996, 

A&A, 312, 599 
Garci'a-Di'az M. T., Henney W. J., Lopez J. A., Doi T., 

2008, Rev. Mex. Astron. Astrof., 44, 181 

Giidel M., Briggs K. R., Montmerle T., Audard M., Rebull 
L., Skinner S. L., 2008, Science, 319, 309 

Heiles C, Crutcher R., 2005, in R. Wielebinski & R. Beck 
ed., Cosmic Magnetic Fields Vol. 664 of Lecture Notes in 
Physics, Berlin Springer Verlag, Magnetic Fields in Dif- 
fuse HI and Molecular Clouds, pp 137 — h 

Henney W. J., Arthur S. J., Garci'a-Di'az M. T., 2005, ApJ, 
627, 813 

Henney W. J., Arthur S. J., Williams R. J. R., Ferland 

G. J., 2005, ApJ, 621, 328 
Herbig G. H., Terndrup D. M., 1986, ApJ, 307, 609 
Hillenbrand L. A., Hartmann L. W., 1998, ApJ, 492, 540 
Hirota T., Bushimata T., Choi Y. K., et al. 2007, PASJ, 

59, 897 

Hogerheijde M. R., Jansen D. J., van Dishoeck E. F., 1995, 

A&A, 294, 792 
Huff E. M., Stahler S. W., 2006, ApJ, 644, 355 
Jansen D. J., Spaans M., Hogerheijde M. R., van Dishoeck 

E. F., 1995, A&A, 303, 541 
Jones B. F., Walker M. F., 1988, AJ, 95, 1755 
Kraus S., Balega Y. Y., Berger J., Hofmann K., Millan- 

Gabet R., Monnier J. D., Ohnaka K., Pedretti E., 

Preibisch T., Schertl D., Schloerb F. P., Traub W. A., 

Weigelt G., 2007, A&A, 466, 649 
Kurucz R. L., 1979, ApJS, 40, 1 

Mellema G., Arthur S. J., Henney W. J., Iliev I. T., Shapiro 

P. R., 2006, ApJ, 647, 397 
Morisset C, Stasihska G., Pena M., 2005, MNRAS, 360, 

499 

O'Dell C. R., Harris J. A., 2010, AJ, 140, 985 

O'Dell C. R., Henney W. J., Abel N. P., Ferland G. J., 

Arthur S. J., 2009, AJ, 137, 367 
Pellegrini E. W., Baldwin J. A., Ferland G. J., Shaw G., 

Heathcote S., 2009, ApJ, 693, 285 
Schilke P., Pineau des Forets G., Walmsley CM., Martfn- 

Pintado J., 2001, A&A, 372, 291 
Schleuning D. A., 1998, ApJ, 493, 811 
Shaw C, Ferland G. J., Henney W. J., Stancil P. C, Abel 

N. P., Pellegrini E. W., Baldwin J. A., van Hoof P. A. M., 

2009, ApJ, 701, 677 



Simon R., Stutzki J., Sternberg A., Winnewisser G., 1997, 

A&A, 327, L9 
Stoerzer H., Hollenbach D., 1998, ApJ, 495, 853 
Stoerzer H., Stutzki J., Sternberg A., 1995, A&A, 296, L9+ 
Tobin J. J., Hartmann L., Furesz G., Mateo M., Megeath 

S. T., 2009, ApJ, 697, 1103 
van der Werf P. P., Stutzki J., Sternberg A., Krabbe A., 

1996, A&A, 313, 633 
Walmsley C. M., Natta A., Oliva E., Testi L., 2000, A&A, 

364, 301 

Wen Z., O'dell C. R., 1995, ApJ, 438, 784 
Young Owl R. C, Meixner M. M., Wolfire M., Tielens 
A. G. G. M., Tauber J., 2000, ApJ, 540, 886 



APPENDIX A: CLOUDY INPUT SCRIPT 

Below we give an example of an input Cloudy script used 
to calculate the grid of photoionization models for the en- 
hanced cosmic rays scenario. 

c 

title - Hydrostatic model of the Orion Bar - 

c 

radius 17.55 18.50 

c 

c Incident continuum 

c 

CMB 

table read "../.. /star_kurucz_39600 . dat " 
Q(H) 49 
brems 6 

luminosity 32.6 range 36.77 to 588.3 Ryd 

c 

c Density and chemical composition 

c 

hden 3.2 

abundances orion 
grains pah 
atom H2 

c 

c Gravity 

c 

constant pressure 
gravity spherical 
gravity external 50 Msun 

c 

c Non-thermal components 

c 

turbulence 2 . km/ s 
magnetic field -5.1 2 
cosmic ray equipartition 

c 

c Stopping conditions 

c 

stop temperature off 
stop eden -4 
set nend 10000 

c 

c Output files 

c 

punch overview "overview.txt" last 
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punch pressure "pressure.txt" last 

punch molecules "molecules.txt" last 

punch lines emissivity "emissivity.txt" last 

S II 6716 

S II 6731 

H 1 6563 

H2 2.121m 

end of lines 

punch grain extinction "extinction.txt" last 

c 

c ... Paranoy@ Ruiz! ; ~D 
c 



